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Abstract: In this talk we shall discuss algorithms and CAD tools for the design and 
analysis of structures for high performance applications using advanced composite 
materials. An extensive mathematical theory for optimal structural (e.g., shape) design 
has been developed over the past thirty years. Aspects of this theory have been used 
in the design of components for hypersonic vehicles and thermal diffusion systems 
based on homogeneous materials. Enhancement of the design methods to include 
optimization of the microstructure of the component is a significant innovation which 
can lead to major enhancements in component performance. Our work is focused on 
the adaptation of existing theories of optimal structural design, e.g., optimal shape 
design, to treat the design of structures using advanced composite materials, e.g., 
fiber reinforced, resin matrix materials. In this talk we shall discuss models and 
algorithms for the design of simple structures from composite materials, focussing on 
a problem in thermal management. We shall also discuss methods for the integration 
of active structural controls into the design process. 
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Problem: Integrated design of structures, their 
materials, and embedded active controls 

Issues: 


1. Shape optimization 

2. Material analysis and design 

3. Actuator design and placement 


Shape Design: Find shape of an object to 
optimize a design criterion and satisfy design 
constraints. 

Abstract Formulation: 
i ncR" the object shape 

• A(u, Q) = 0 defines u(x ) € R m , * G ^ 

• Given /(u,Q) a real-valued function 
Optimal Shape Design Problem: 


Two essential problems: 


1. Select the topology for the structure (cylin- 
drical, rectangular, etc.); and 


2. Within the designated topology find the 
best shape. 


Remarks: 


• The first problem is very difficult; e.g., in- 
troduction of internal holes in a structure 
to reduce the weight without violating de- 
sign constraints 

• The second problem (initial and final topolo- 
gies are the same) can usually be treated 
by gradient methods. 


Example: Optimal Compliance Design of an 
Elastic Structure 

Problem: Design an elastic structure contain- 
ing a large number of "cells" in a continuous 
array; e.g., fiber reinforced structure. 

Remark: If the array is locally periodic, the 
macroscopic moduli may be computed using 
homogenization theory. 

Design Parameters: dimensions (a, 6) and ori- 
entation 9 of the microscopic elements. 
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Design Parameters 



Design Algorithm: 


1. Use homogenization to compute the local 
effective elasticity tensor E H (x, (a, b, c, 0)). 

2. Compute gradient of performance func- 
tion 


3. Steepest descent on design parameters 



Heterogeneous Structure 
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Effective Parameter Model: 

“Effective” elasticity tensor E H (x ) may be 
computed using homogenization theory. 


1. Solve cell problem 
2 


ijmn - 


E .. 9 Xm° dt'i, _ 

=, JY E ' jmn ~diM dy 3 ~ 


£ J 


_ dv{ • 
_,JY Ei i mn dy? Y 


t,y t m,n=l ’ 

for the characteristic deformation Xm^ 


2. Compute the homogenized elasticity ten- 
sor 
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Displacement Model: Assume smooth varia- 
tion across structure. 


t?H, (a, b, 9) 
h ijkl 


(x) 


2 


E 

IJKL - 1 




Rjj [0 C * )] RkK [0 ( x )] R IL (* )] 


*[*(x)] 


cos 9 -sin0 
sin0 cos0 


is the (local) rotation matrix. 


Macroscopic Behavior: (For each © — (a, b,9)) 


Ll 


ijkl= 1 


ijki 


(x) 


dxi dxj 


= ti n ^‘ dn+ ^k Tmdr 

fi,i = 1,2 are the applied body forces in £2, 
r .j = 1,2 are the tractions applied on the 
boundary f T C r = 9Q. 



Optimal Compliance Design 


min V / 

(®=(a,b,0)} Jn l > kl dx t dxj 
subject to the constraint (*) and 


f (1 - a(x)b(x))dQ. < |f2p| 

Cl F = the maximum volume fraction allocat- 
ed to the reinforcing material. 


Using a penalty method, the optimization prob- 
lem is approximated by 

max min n £ (i/) 

© v£V 

where n £ <» is the total potential energy 
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Introduce the Lagrangian 


£ = IT - A( / (1 — ab)dQ — ftj?) 

JO. 

where A < 0 is a Lagrange multiplier. 


Taking the variation of £ with respect to u 
and the design variables © = (a, 6,0) gives 
the optimality conditions: 


biju = i 


1 

2 


du k dui 
da dxi dx j 


f A b 


Sadn > 0 


V6a = a* — a, 0 < a* < 1 E Q 
* * * 


Optimality conditions : 


a r= min < max 


ft = min < max 


0 , d — Pa 


o ,ft - Pi, 


/ ATT* 1 *® 

^2 da dxidxj 

1 dukdui 


+ Aft 


,1 


Aa 


^2 db dx t dxj 

A = min { 0, A - p\ {^J (1 - ab)dSl - j 


1 ° E ijkt du k dui 


= 0 


2 dO dx t dxj 
for arbitrary positive numbers p a ,Pb>P\ 



Example: Optimal design of a thermal diffus- 
er using composites 

Problem: Select the shape Q of the diffus- 
er and the parameterization © of S2f (vol- 
ume fraction, orientation, packing, etc.) of 
the material infrastructure to minimize the 
weight of the diffuser and meet operational 
objectives. 

(i) the maximum temperature at the payload 
- diffuser interface must not exceed T m ; 

(ii) no part of the diffuser can be thinner than 
some constant d; and 

(iii) for convective cooling, the flux on the 
transmission interface of the diffuser must be 
below q m . 


Composite Material Thermal Diffuser 
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Model: 


Assume: Conductivity of the fiber material is 
k, conductivity of the matrix material is K 
(could be anisotropic) 

Conductivity tensor: {a9(x,y,z),i,j = 1,2,3} 


d 

dx± 


‘§ ( *>O =0 


with the boundary conditions 

dT@ =qin on r i; I^ = 0on r 2 


dr ) a @ 


dv a @ 


dT 


.© 


dv„@ 


*(T @ ) P = g, on Z 3 


d/dr) a& is the conormal derivative at the sur- 


face. 


Design Parameters: ©, the fiber orientation 
and packing, and R(z), 0 < z < L the curve 
defining the shape of the boundary, and L. 

Performance Index: Mass of diffuser 

n(©, £(•),£) = *• [ if p @ (x,y,z)R(z) 2 dxdydz 
JO J Jc(z) 

p®(x,y,z ) is the mass density in a cross sec- 
tion C(z ). 

The optimal design problem is 
min n(©, R(-)> L) 
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Homogenization: (local) 


Effective Conductivity: 


a ij = M(a {j ) - M 


Macroscopic Behavior: 


^Au = - au J^- = f 


tJ dxjdxj 


‘Corrector:" 


_ d ( , . a j\ _ da ij(y ) 

dyi { ' }(y) d yj x ) d yi 


Optimal Design: 


Adjoin constraints to performance function: 
n „(©, R(-), l)(v) = n(©, l ) 

M > 0 is a small parameter. 
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Define the Lagrangian: 


C(v,&,R(-),L,A) = 

— A[ / ppdxdydz - Mp ] 

JQ]? 

Lagrange multiplier A < 0. 

Optimality Criteria : 

Obtained from variation of C with respect to 
state u and design variables (©, R(-), L). 


